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ABSTRACT 


This  report  illustrates  a  proposed  definition  of  surface  density 
of  contamination  for  the  deposit  of  fractionated  fallout.  The  illus¬ 
trations  employed  give  estimates  of  the  dependence  of  exposure-dose 
rate  on  the  degree  of  radionuclide  fractionation  and  of  the  sensitivity 
of  this  dependence  to  independent  fission  yield  input  data,  total  yield 
of  the  device,  and  time  of  exposure. 
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SUMMARY 


Problem 


The  effects  of  radionuclide  fractionation  severely  complicate  the 
prediction  of  many  properties  of  nuclear  bomb  debris,  including  the 
definition  of  surface  density  of  contamination  and  the  dependence  of 
exposure -dose  rate  on  the  degree  of  fractionation. 

Findings 

A  semi -empirical  model  can  be  used  to  illustrate  a  nev  definition 
of  surface  density  of  contamination  for  fractionated  fallout.  The 
model  lo  recommended  for  rule-of-thumb  estimates  of  the  effect  of 
fractionation  on  exposure-dose  rate,  und  us  a  stop-gup  until  either 
better  models  or  more  extensive  information  becomes  available. 
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PREFACE 


This  series  of  reports  presents  and  discusses  the  effects  of 
radionuclide  fractionation  in  nuclear  bomb  debris.  Bart  I  (Reference 
1)  defined  fractionation  as  "any  alteration  of  radionuclide  composition 
occurring  between  the  time  of  detonation  and  the  time  of  radiochemical 
analysis  which  causes  the  debris  sample  to  be  nonrepresentative  of  the 
detonation  products  taken  as  a  whole."  It  showed  how  the  radionuclide 
compositions  of  fractionated  camples  could  be  correlated  empirically 
by  logarithmic  relations.  Part  II  (Reference  2)  UBed  these  relations 
as  the  basis  of  a  technical  discussion  of  contamination  density  as 
applied  to  fractionated  nuclear  debris.  Part  III  (Reference  3)  presen¬ 
ted  a  theoretical  foundation  for  the  observed  logarithmic  correlations 
of  rart  I.  It  used  this  as  a  simplified  meanG  of  estimating  fraction¬ 
ation  as  a  function  of  particle  size  and  Lhe  partition  of  product  radio¬ 
nuclides  among  local,  intermediate,  and  worldwide  fallout. 

The  present  report  extends  the  calculations  of  Ihrt  III  to  show  how 
fractionation-correlation  parameters  can  be  used  to  estimate  the  expo¬ 
sure-done  rate  from  nuclear  debris  with  various  degrees  of  fractionation. 
It  serves  to  illustrate  the  proposals  made  in  Bart  II, 
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INTRODUCTION 


Local  fallout  from  a  land-surface  nuclear  detonation  consists  of 
the  larger  radioactive  particles  formed  in  the  explosion.  Because  of 
radionuclide  fractionation,  it  is  depleted  in  volatilely  behaving  mass 
chains  (such  as  the  mass -89  chain)  relative  to  refractorily  behaving 
mass  chains  (such  as  the  mass-95  chain).  Worldwide  fallout,  on  the 
other  hand,  consists  of  the  smaller  particles  of  debris  and  is  rela¬ 
tively  enriched  in  volatilely  behaving  mass  chains.  The  effect  of 
fractionation  on  exposure-dose  rate  from  fallout  deposits  is  to  pro¬ 
duce  a  wide  variation  in  the  ratio  of  dose-rate  to  f is si on -product 
abundance,  for  any  given  mass  chain,  or  of  one  fisBi on-product  abund¬ 
ance  tu  another,  for  certain  pairs  of  mass  chains. 

The  radiochemical  composition  of  fallout  will  depart  from  the 
representative  composition  of  debris  by  a  quantity  which  we  have  called 
the  degree  of  fractionation  and  have  defined  quantitatively  as  the 
bane-10  logarithm  of  the  fractionation  ratio: 


*  a  log10  r89,95 


The  fractionation  ratio  (rgg  95)  is  the  ratio  of  the  number  of  fissions 
(from  the  device  in  question)  required  to  produce  the  quantity  of  mass- 
89  cliain  found  at  a  given  location,  to  the  number  required  to  produce 
the  quantity  of  mnss-95  chain  at  that  location.  The  relation  of  our 
notation  to  that  of  other  authors  is  given  in  Appendix  A. 

Two  theoretical  methods  exist  for  estimating  $  in  local  fall¬ 
out,  3, f*,5,o  empirical  relationships  are  available  for  estimating 

the  relative  amounts  of  other  mass  chains  as  a  function  of  $.1  By  com¬ 
bining  either  of  these  theoretical  methods  with  the  empirical  relations, 
one  obtains  a  semierapirical  means  of  estimating  the  radionuclide  compo¬ 
sition,  and  hence  the  resulting  exposure-dose  rate,  as  a  function  of 
particle  size.  A  fallout-transport  model  can  then  be  used  to  calculate 
these  same  quantities  at  any  point  in  the  fallout  pattern. 
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The  purpose  of  this  report  is  threefold:  (l)  to  illustrate  the 
empirical  estimation  of  radionuclide  composition  as  a  function  of  4; 

(2)  to  estimate  the  resulting  exposure-dose  rate  as  a  function  of  *, 

{3;  to  obtain  a  preliminary  estimate  of  the  sensitivity  of  this  relation¬ 
ship  to  several  kinds  of  input  data;  and  (h)  to  illustrate  the  applica¬ 
tion  of  a  new  definition  of  surface  density  of  contamination. 2  In  pre¬ 
senting  this  material,  familiarity  with  the  previous  reports  in  this 
series  is  assumed. 

The  procedure  for  estimating  the  radiochemical  composition  of  a 
fission-product  mixture  as  a  function  of  4  is  as  follows.  We  consider 
only  those  nuclides  which,  at  some  time,  contribute  significantly  to 
the  dose  rate.  These  nuclides  and  their  isobars  form  the  significant 
mass  chains.  Next,  the  fission-product  elements  are  divided  into  re¬ 
fractory  and  volatile  groups,  according  to  their  vapor  pressures  at  a 
chosen  temperature.  The  temperature  chosen  in  this  report  corresponds 
roughly  to  the  solidification  temperature  of  molten  silicate  soil.  For 
any  nuclear  detonation,  the  total  yield  of  the  device  determines  the 
solidification  time.*  In  any  chain,  the  fraction  of  atoms  that  exists 
in  the  form  of  refractory  elements  at  the  time  of  solidification  (Fr) 
is  a  measure  of  tlv  refractory  nature  of  Mint  nnnn  rhnln  taken  an  a  whole. 
According  to  the  magnitude  of  the  fraction,  the  chains  are  grouped  into 
three  categories:  volatile,  refractory,  or  mixed.  The  specification  of 
o  ( I .e. ,  the  degree  of  fractionation)  will  be  shown  to  fix  the  value  of 
the  fractionation  coefficient  (r^  gr()  for  ench  chain  k  and,  therefore, 
the  radiochemical  composition. 

Part  II  of  this  series2  showed  tliat  one  could  estimate  the  effect 
of  fractionation  on  exposure -dose  rate  by  summing  the  products  of  the 
f i-actionation  coefficients  and  the  unfractionated  contributions  for  the 
individual  chains.  In  so  doing,  it  was  convenient  to  divide  the  chains 
into  groups,  according  to  the  volatility  shown  by  their  behavior.  This 
procedure  will  be  carried  out  In  detail  in  the  present  report.  Thus, 
these  calculations  will  illustrate  the  definition  of  surface  density 
of  contamination  proposed  in  Part  II. 


♦The  energy  release  from  a  nuclear  device  is  given  sometimes  as  the 
yield  from  fission  reactions  and  sometimes  as  the  total  (fireball) 
yield.  The  latter  includes  the  former  plus  the  energy  released  by 
fusion  reactions.  The  yield  from  fission  reactions  plays  the  domin¬ 
ant  role  in  producing  radioactivity.  The  total  or  fireball  yield, 
together  with  the  nature  and  location  of  environmental  material, 
determines  the  condensation  time  and  solidification  time  for  the  debris. 
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PRELIMINARY  CONSIDERATIONS 


Selection  of  Significant  Chains 

Although  some  90  mass  chains  are  produced  in  the  thermal  neutron 
fission  of  u235,  only  about  half  of  these  will,  at  some  time  or  other, 
make  a  significant  contribution  to  the  exposure  dose  rate.  Therefore, 
calculations  can  be  considerably  shortened  by  considering  only  signi¬ 
ficant  mass  chains.  In  selecting  those  chains,  and  in  all  subsequent 
dose-rate  calculations,  the  contributions  of  individual  fission  products, 
such  as  those  compiled  by  Miller  and  Loeb  for  the  products  of  thermal- 
neutron  fission  of  U235}  vill  be  used.*?  For  any  given  time,  the  nuc¬ 
lides  contributing  less  than  0.3  $  to  the  total  dose  rate  at  that  time 
will  be  considered  insignificant.  A  significant  nuclide  and  its  isobars 
form  a  significant  maos  chain.  These  47  significant  chains  will  be 
listed  in  Table  1,  together  with  the  time  interval  in  which  they  are 
significant.  The  earliest  time  considered  will  be  1+5.8  m. 

Criterion  of  Volatility 

The  empirical  .  •ItiLloiishipn  derived  in  Reference  1  depend  upon  a 
quantity,  Fr,  equal  to  the  fraction  of  n  particular  chain  that  is  re¬ 
fractory  at  the  time  of  condensation  in  the  fireball.  Miller's  thermo¬ 
dynamic  model' 1 for  land -surface -burst  fallout  adopts  an  idealized 
soil  that  melts  at  l400°C.  Above  this  temperature,  Miller's  first  stage 
of  condensation  pertains,  and  the  product  nuclides  distribute  themselves 
uniformly  throughout  the  volume  of  molten  particles.  Below  this  tempera¬ 
ture  the  particles  are  frozen  and  produce  nuclides  can  only  deposit  on 
the  surfaces.  Miller  calls  this  the  second  stage  of  condensation.  In 
order  to  maintain  consistency  between  the  two  approaches,  we  will  also 
adopt  this  approximation  and  call  nil  elements  volatile  whose  predomin¬ 
ant  species  liave  a  normal  boiling  point  less  than  l400°C.  These  ele¬ 
ments  are  As,  Sc,,  Br,  Kr,  Rb,  Mo,  Tc,  Te,  I,  Xe,  and  Cs.  All  other 
elements  will  be  called  refractory. 

Calculation  of  Fr 


To  calculate  F^,  we  now  need  an  equation  to  estimate  the  time  at 
which  l400°C  Is  reached  by  the  cooling  fireball.  From  Table  3.8  of 
Reference  6,  the  equation 

t  (see)  =  1.88  w0'363 


can  be  inferred,  where  V  is  the  total  yield  in  kt. 
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TABLE  1 


Fj  Values  for  Significant  Chains 


Ctaln  6-»«c  SoUdlTlcbtlcn  4l-»«o  Solidification  Tit  of  Sl*nlflcano«» 

Vo.  RwS  8Sn&5E  RSamE  01«ndonln  froa  fo 


$ 

0.01 

0.04 

0.00 

0.00 

65.8  a 

1.12  b 

is 

0.00 

0.01 

0.00 

0.00 

65.8  a 

3.52  h 

87 

0.00 

0.00 

0.00 

0.00 

2.60  h 

23.8  h 

0.00 

0.00 

0.00 

0.00 

65.8  a 

n.i  h 

88 

0.00 

0.00 

0.00 

0.00 

65.8  a 

23.8  b 

8? 

0.00 

0.00 

0.00 

0.00 

6s. 8  a 

2.6o  h 

90 

0.01 

0.01 

0.08 

0.08 

2.63  y 

55-3  y 

91 

0.06 

0.07 

0.16 

0.17 

65,8  a 

3.12  d 

92 

0.28 

0.30 

0.82 

0.83 

65.8  a 

1.6s  d 

93 

O.fiO 

0.62 

0.99 

0.99 

65.0  a 

3.12  d 

9* 

0.88 

0.87 

0.99 

0.99 

65.8  a 

e.4o  d 

95 

97 

0.97 

0.97 

1.00 

1.00 

16.2  h 

2.60  y 

1.00 

1.00 

1.00 

1.00 

65. 8  a 

6.70  d 

98 

1.00 

1.00 

1.00 

1.00 

65.8  a 

3-52  »> 

99 

0.99 

0.97 

0.86 

0.83 

7.56  h 

21.1  d 

101 

0.O0 

o.do 

0.09 

0.07 

65.8  a 

2.40  h 

102 

0.61 

0.39 

0.02 

0.01 

63.8  a 

1.&  h 

103 

0.3** 

0.20 

0.17 

0.17 

26.2  1, 

1.20  y 

105 

0.12 

0.21 

0.75 

0.80 

65.8  m 

4.57  d 

106 

0.20 

0.36 

0.93 

0.95 

66.4  d 

0.18  y 

107 

0.36 

0.56 

0.99 

0.99 

65.8  » 

45.8  a 

118 

1.00 

1.00 

1.00 

1.00 

1.20  y 

55.3  y 

125 

i.cc 

1.00 

1.00 

1.00 

1.20  y 

12.0  y 

126 

1.00 

1.00 

1.00 

1.00 

2.40  h 

2.13  d 

127 

1.00 

1.00 

1.00 

1.00 

2.13  a 

14.4  d 

128 

1.00 

1.00 

1.00 

1.00 

1.66  h 

4.57  d 

129 

0.97 

0.98 

0.97 

0.98 

65.8  a 

23.8  h 

131 

0.74 

0.86 

0.73 

0.86 

65.8  a 

66.4  d 

132 

0.48 

0.70 

0.60 

0.58 

1.12  h 

30.9  d 

133 

0.25 

0.53 

0.23 

0.50 

65.8  n 

30.9  d 

13<> 

0.09 

0.26 

0.01 

0.02 

65.8  a 

11.1  b 

135 

0.00 

0.10 

0.00 

0.00 

65.8  a 

3.12  d 

137 

0.00 

0.00 

0.00 

0,00 

66.4  d 

250  y 

138 

0.00 

0.00 

0.00 

0.00 

45.8  a 

5.26  h 

139 

0.00 

0.02 

0.06 

0.06 

45.8  a 

1.11  h 

Ho 

0.06 

0.11 

0.28 

0.32 

4.56  b 

H3  d 

Hi 

0.38 

0.k7 

0.96 

0.96 

63.8  a 

301  d 

H2 

0.71 

0.78 

1.00 

1,00 

45.8  a 

1.11  b 

H3 

0.94 

0.96 

1.00 

1.00 

65.8  a 

H.4  d 

H4 

0.98 

0.99 

1.00 

1.00 

H.4  d 

8.18  y 

H5 

1.00 

1.00 

1.00 

1.00 

65.8  a 

2. 13  d 

H6 

1.00 

1.00 

1.00 

1.00 

65.S  a 

3.52  b 

H? 

1.00 

1.00 

1.00 

1.00 

23.8  b 

97.3  d 

H9 

1.00 

1.00 

1.00 

1,00 

65.8  a 

H.4  d 

151 

1.00 

1.00 

1.00 

1.00 

7.56  h 

4.57  d 

152 

1,00 

1.00 

1.00 

1.00 

45-6  a 

3.52  h 

155 

1.00 

1.00 

1.00 

1.00 

1.78  y 

3-80  y 
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Bo lies  and  Ballou®  tabulate,  for  U83 5 -thermal  neutron  fiBaion,  the 
number  of  atoms  of  each  fission-product  masB  chain  present  in  various 
elemental  forms  at  these  times.  These  elemental  distributions  are 
calculated  both  according  to  Present's  minimum  kinetic  energy  (MKE) 
theory  of  charge  distribution, 9  and  according  to  the  equal  charge  dis¬ 
placement  (ECD)  theory  of  Glendenin,  Coryell  and  Edwards.10  Although 
the  latter  treatment  is  more  widely  preferred,  we  will  cany  out  cal¬ 
culations  on  both  bases  in  order  to  illustrate  the  sensitivity  of  the 
calculations  to  the  fractional  chain  yields.* 

A  sample  calculation  of  Fr  for  the  mass-92  chain,  according  to 
Present's  theory,  for  a  solidification  time  of  hi  sec, is  given  below. 
Values  of  Fr  for  the  significant  mass  chains  are  given  in  Table  1  for 
both  theories  of  independent  yield  (nuclear  charge  distribution)  for 
solidification  times  of  6  and  hi  sec.  According  to  the  equation  above. 


these  correspond  to  total  yields 

Mass-92  Decay  Chain  Half-Life 

of  about  25  kt  and 

Atans  per  10^ 
Fissions  at 
hi  sec.  (Ref. 8) 

5  Mt,  respectively 

Refractory  Atoms 
per  10’!  FitiBiuns 

Br 

1.5  sa 

0 

0 

i 

Kr 

3  c 

0 

0 

i 

Rb 

| 

17  oa 

102 

ob 

4 

Sr 

2.7  h 

h77 

h77 

l 

Y 

3.5  h 

1 

1 

Zr 

stable 

0 

0 

Total 

580 

478 

fr  =  478/500  =  0.82 

a.  Estimated  half-life. 

b.  For  silicate  soil,  Rb  would  probably  behave  refractorily  because 
it  forms  refractory  silicate  compounds. 


*The  fractional  chain  yield  of  a  fission-product  radionuclide  is  the 
ratio  of  the  independent  (primary)  yield  of  that  nuclide  to  the  total 
chain  yield  for  that  nuclide's  mass  chain.  The  concept  applies  only 
at  the  instant  of  fission. 
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Volatile,  Mixed,  and  Refractory  Chains 


The  significant  mass  chains  can  now  be  divided  into  three  groups 
(volatile,  mixed,  or  refractory)  according  to  the  Fp  value  for  that 
chain.  Chains  with  O.98  <  Fp  <  1.00  are  called  refractory,  chains  with 
0.02  <  Fp  <  O.98  are  called  mixed,  and  chains  with  0.00  Fp  ^  0.02 
arc  called  volatile. 

For  unfractionated  debris,  the  total  exposure -dose  rate  und  the 
contributions  of  the  individual  groups  are  shown  as  functions  of  time 
in  Figs.  1  through  4  for  the  two  yield  theories  and  the  two  condensa¬ 
tion  times.  The  total  curve  is,  of  course,  the  same  in  all  four  fig¬ 
ures  . 


Hie  percentage  contributions  of  each  group  to  the  total  at  1.12 
and  23.8  hr  are  shown  in  Table  2.  As  can  be  seen,  98  $  or  more  of  the 
dose  rate  la  accounted  for  by  the  significant  chains. 

TABLE  2 

Contribution  of  Groups  to  Unfractionated  pose  Rate  ($) 


Group 

1 .12 -hr 

close  rate 

23.8-ln 

duoo  rule 

Present 

Glendenin 

Present 

Glendenin 

6- sec 

Solidification 

•Time 

Volatile 

33-4 

30.9 

16.8 

0.4 

Mixed 

47.0 

49.2 

91.0 

67.O 

Ref  ractory 

18.5 

18.0 

30.4 

29.9 

Neglected 

1.1 

1.1 

1.8 

1.9 

4l-sec  Solidification 

Time 

Volatile 

48.4 

48.4 

16.8 

16.8 

Mixed 

21.2 

20.9 

38.6 

38.0 

Ref rectory 

29.2 

29.5 

42.7 

43.3 

Neglected 

1.2 

1.2 

1.9 

1-9 
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EXPOSURE  OOSE  RATE  (nr/hf) 


TIME  (OAYSI 


Fig.  1  Exposure-dose  Rates  From  Volatile,  Refractory,  Mixed,  and  Total 
Fission-Product  Mass  Chains  as  Functions  of  Time.  Calculated 
for  3  ft  above  an  infinite,  smooth,  impenetrable  plane,  uniformly 
contaminated  with  the  products  of  104  thermal-neutron-induced 
fissions  of  U235  per  sq..  ft.  (6-sec  solidification.  Present's 
Theory . ) 
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6-SEC  SOLIDIFICATION 
(Glendenin's  Theory) 


TIME  (DAYS) 


Exposure -dose  Rates  From  Volatile,  Refractory,  Mixed,  and  Total 
Fission -Product  Mass  Chains  as  Functions  of  Time.  Calculated 
for  3  ft  above  an  infinite,  smooth,  impenetrable  plane,  uniformly 
contaminated  with  the  products  of  lcr  thermal-neutron-induced 
fissions  of  U^35  per  gq.  ft,  (6-sec  solidification,  Glendenin's 
Theory.) 
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0  01  01  10  10  100  1000 

TIME  (DAYS) 


Fig.  3  Exposure-dose  Rates  From  Volatile,  Refractory,  Mixed,  and  Total 
Fission-Product  Mass  Chains  as  Functions  of  Time.  Calculated 
for  3  ft  above  an  infinite,  smooth,  impenetrable  plane,  uniformly 
contaminated  with  the  products  of  10^  thermal-neutron-induced 
fissions  of  U235  per  sq,.  ft.  (4l-sec  solidification,  Present's 
Theory. ) 
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TIME  (DAYS) 


Fig.  4  Exposure -dose  Rates  From  Volatile,  Refractory,  Mixed,  and  Total 
Fission -Product  Mass  Chains  as  Functions  of  Time.  Calculated 
for  3  ft  above  an  infinite,  smooth,  impenetrable  plane,  uniformly 
contaminated  with  the  products  of  10*  thermal-neutron-induced 
fissions  of  U235  per  sq.  ft.  (kl-sec  solidification,  Olendenin's 
Theory . ) 
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RELATION  TO  THE  NEW  DEFINITION  OF  CONTAMINATION  IEVEL 


i 


Part  II  of  this  series'"  proposed  a  new  definition  of  surface  den¬ 
sity  of  contamination  which  was  applicable  to  fractionated  debris.  The 
new  definition  was  summarized  by  Eq.  18  of  that  report  as 

1"bk 

»<‘>  "  “95  [s  W*>  *  1  **  E 

quasi-refractory  quasi -volatile  induced 

chains  chains  chains 


where  Onr  is  the  contamination  surface  density  in  mass -95  chain  equi¬ 
valent  fissions  per  unit  area.  is  the  average  number  of  fission- 
product  or  induced  atoms  of  mass  k  produced  initially  per  fission,  and 
Djc(t)  is  the  dose  rate  contribution  of  the  maBB-k  chain  per  atom  of 
mass  k  per  unit  area  at  the  time  t.  The  quantity  b^  is  a  correlation 
parameter  which  is  dicussed  in  Ref.  (3),  and  which  has  been  observed  to 
be  approximately  equal  to  Fp1/'^  in  debris  from  high-yield  surface 
bal’d  Ls .  -  >  * 

For  illustration  and  interim  predictions,  the  "air  ionization  rates" 
(actually  the  exposure-dose  rates)  of  individual  fission-product  radio¬ 
nuclides  tabulated  by  Miller  and  Loeb  provide  a  convenient  source  of 
data.  These  authors  tabulate  in  units  of  nr/hr,  the  "air  ionization 
rate"  three  feet  above  an  infinite  plane,  each  square  foot  of  which  is 
uniformly  contaminated  with  the  products  of  l(r*  thermal-neutron-induced 
fissions  of  Us35,  in  their  Ihble  A-l,  they  list  the  Individual  contri¬ 
bution  of  each  fission -product  radionuclide  at  various  convenient  times. 
In  the  notation  of  Ref.  2,  (cf.  Eq-  6)  each  entry  would  be  written 


j 


with  the  contamination  surface  density  (o^)  in  units  of  10^  fissions 
per  sq  ft,  the  total  chain  yield  for  mass  k  (Y^)  being  the  value  for 
thermal -neutron  fission  of  U^35,  and  the  product  cf  the  dose-rate  con¬ 
version  factor  (dj)  the  chain  fraction  and  the  decay  constant 

(kj)  being  in  units  of  nr  per  hr  per  fission.  Since  the  tabulated  data 

♦This  relation  neglects  the  departure  from  zero  of  the  term  ak  of  Refa. 
2  and  3.  When  dealing  with  samples  of  fallout,  rather  than  particles 
of  a  single  size,  this  is  Justifiable. 


11 


is  for  unfractionated  fission  products,  each  radionuclide  teis  the  same 
value  of  0^.  Thus ,  Dj  gives  the  dose  rate  from  each  radionuclide  under 
the  conditions  stated,  and  if  the  values  are  sunned  for  all  nuclides 
of  a  given  mass  k,  one  obtains  D^,  the  dose  rate  from  the  mass-k  chain 
at  the  time  t.  In  Ref.  2,  it  was  found  convenient  to  abbreviate  the 
sum  EjdjGjj^t)^  by  D^(t)  and  write  (Ref.  2,  Eq.  ^) 


Dk  =  °kYk\(t) 


To  illustrate  the  application  of  Eq.  18  with  Miller  and  Loeb's 
data,  we  will  consider  first  the  quasi-refractory  term.  The  value  of 
this  term  is  estimated  by  (l)  choosing  from  Table  1  the  mass  chains 
with  1.00  Fr  ^0.98  for  the  appropriate  conditions,  (2)  summing  all 
the  data  tabulated  by  Miller  and  Loeb  for  these  chains  at  the  time  of 
interest  t. 

In  calculating  the  value  of  the  qua si -volatile  term,  it  has  been 
found  convenient  and  instructive  to  consider  the  mass  chains  to  fall 
into  two  classes:  n  purely  volatile  class  with  0.02  ^  Fr  ^.0.00  and  a 
mixed  class  with  O.98  >  Fr  >  0.02,  For  the  purely  volatile  class,  we 
(i)  choose  the  appropriate  mass  chains  r.r  above,  (2)  sum  their  contri¬ 
butions  for  the  time  of  interest,  und  (3)  multiply  the  sum  by  the  frac¬ 
tionation  ratio  we  wish  to  apply.  The  contribution  for  the  mixed  chains 
is  calculated  similarly,  except  that  in  each  case  the  fractionation  1  , 
ratio  muot  be  raised  to  the  appropriate  power  l-bk,  here  taken  as  1  -  F^  '  , 
for  the  chain  in  question. 

The  induced  chains  are  not  considered  in  this  report.  The  value 
of  0513  is  taken  as  10*  fissions  per  square  foot. 
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CALCULATIONS  AND  RESULTS 


Exposure -Dose  Rates 

Each  term  in  the  summation  of  D(t)  can  be  considered  as  the  pro¬ 
duct  of  the  contribution  from  unfractionated  fission  products 
(°95  and  a  fractionation  correction  factor  (rgg  05)^“^.  Fig¬ 

ure  5  illustrates  the  variation  of  this  correction  factor  with  frac¬ 
tionation  ratio  (^9,95)  as  computed  for  several  values  of  Fr. 

As  a  specific  example,  consider  the  mass-132  chain  contribution 
at  1.12  h  from  a  burst  with  a  4l-sec  solid  1ft cation  time.  From  Itablc  1, 
the  value  of  Fr  according  to  Present's  theory  is  0.40.  Fran  Fig.  5» 
the  correction  factors  for  ^9^95  values  of  l/3,  l/lO,  l/30  and  l/lOO 
would  be  O.667,  0.427,  0.264  and  O.lOC  respectively.  The  unfraction¬ 
ated  contribution  of  this  chain  to  the  total  dose  rate  is  calculated 
from  Ref.  7  as  0.01010  nr/hr/lO^  fissions/ft''-.  Contributions  which 
the  chain  would  in  fact  provide, for  the  fractionation  ratios  listed, 
are  therefore  0.00674,  0.00431,  0.00287  nnd  0.00184  nr/hr/lO+  fissions/ 
ft^,  respectively.  Similar  calculations  of  contributions  for  all  the 
other  mixed  chains  and  summation  of  these  contributions  give  values  of 
0.846,  0.554,  0.434  and  0.347  ni'/hr/lO1*  fissions/ft2  respectively,  at 
these  fractionation  ratios.  Thbles  3  through  6  liot  the  contributions 
of  each  group  for  each  set  of  conditions  considered,  together  with 
their  totals.  They  are  given  both  in  units  of  nr/hr  for  10^  mass-95 
chain  equivalent  fiseions^sq  ft  and  in  kr/hr  for  1.45  x  1023  mass-95 
chain  equivalent  fissions/sq  mi.  To  illustrate  the  functional  behavior 
of  the  contributions,  Figs.  6  through  9  shew  the  variations  of  compon¬ 
ent  contributions  and  total  dose  rates  with  the  fractionation  ratio,  so 
that  the  effects  of  solidification  time,  independent-yield  theory  and 
exposure  time  can  be  independently  observed.  The  figures  are  drawn 
for  a  constant  number  of  masG-95  chain  fissions.  On  this  basis  the 
contribution  of  the  refractory  group  (bj  =  l)  is  constant  and  the  contri¬ 
bution  of  the  volatile  group  (bfc=0)  is  directly  proportional 

*The  number  of  mass-95  chain  equivalent  fissions  (fo^)  is  equal  to  the 
number  of  device  fiGsions  which  produced  the  quantity  of  mass-95  chain 
observed. 
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Fractionation  Correction  Factors  as  Functions  of  F„  and  rar. 
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Difference  due  to  rounding -off  errors. 
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Differences  due  to  rounding  off  errors 


Fig.  6  Relative  Exposure-Dose  Rate  From  Volatile-Depleted  Fission 

Products  as  a  Function  of  Fractionation  Ratio.  Same  geometry 
as  E'lg.  1.  (1.12  h,  6-sec,  solidification,  Glendenin's  Theory) 


f  R'J.S'j 


Fig.  7  Relative  Exposure-Dose  Rate  Fran  Volatile -Depleted  FisBlon 

Products  as  a  Function  of  Fractionation  Ratio.  Same  geometry 
as  Fig.  1.  (1,12  h,  4l-sec,  solidification,  Glendenin’s  Theory) 
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Fig.  8  Relative  Exposure-Dose  Rate  From  Volatile-Depleted  Fission 

Products  as  a  Function  of  Fractionation  Ratio.  Rams  geometry 
as  Fig.  1.  (23.8  h,  4l-sec.  solidification,  Glendenin's  Theory) 


Fig.  9  Relative  Exposure-Dose  Rate  From  Volatile -Depleted  Fission 

Products  as  a  Function  of  Fractionation  Ratio.  Same  geometry 
as  Fig.  1.  (23.8  h,  tl-sec.  solidification,  Present's  Theory) 
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to  the  fractionation  ratio.  In  order  to  include  the  contribution  of 

induced  activities,  one  need  only  increase  the  contribution  of  the  i 

refractorily  behaving  group  by  the  amount  appropriate  to  the  device  '* 

and  environmental  conditions  of  interest. 

Estimated  Decay  and  the  Way-Wigner  Rule 

By  assuming  that,  at  reasonably  early  times  unfractionated  mixed 
fisaion  products  behave  as  a  statistical  assembly  of  emitters  with 
decay  constants  linearly  related  to  the  fifth  power  of  the  disintegra¬ 
tion  energy,  Way  and  Wigner11  derived  a  t-1*2  dependence  for  the  gross 
beta  decay.  The  rule  has  been  found  useful  in  estimating  the  exposure- 
dose  rate  decay  of  mixed  fission  products, but  has  become  the  subject 
of  a  popular  fallacy,  namely,  that  departures  from  the  rule  indicate 
fractionation,  while  correspondence  to  the  rule  indicates  representa- 
tivity.  To  illustrate  the  effect  of  fractionation  on  the  decay  of 
exposure  dose  rate,  Table  7,  lists  the  ratio  of  the  1.12-hr  exposure- 
dose  rate  to  the  23.8-hr  exposure-dose  rate  as  calculated  for  the 
various  conditions  listed  here  and  from  the  t"-*-*2  rule.  The  calcula¬ 
tions  show  that,  rather  than  produce  a  departure  from  the  t“l*2  rule, 
it  is  possible  that,  in  the  absence  of  induced  activities,  fractiona¬ 
tion  could  promote  correspondence  to  the  rule.  The  main  point  here  is 
to  beware  of  superficial  generalizations  based  on  ouch  a  complicated 
property  of  nuclear  debris  as  gross  decay  rate. 


TABLE  7 

Ratio  of  1.12-hr  Expoaure-Dose  Rate  to  23-8-hr  Exposure  Dose 
Rate  for  Various  Conditions 


r89,95 

6-sec  Solidification 
Present  Glendenin 

4l-sec  Solidification 
Present  Glendenin 

1/3 

46 

46 

43 

42 

1/10 

4l 

4l 

38 

38 

1/30 

40 

4o 

36 

36 

1/100 

40 

39 

37 

36 

0 

32 

34 

37 

36 

Unfractionated  value  (?Qg  g^=l) : 

53 

Way-Wigner  value  (1.12/23.8)'!’2: 

39 
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DISCUSSION 


Although  the  calculations  described  here  are  primarily  intended 
to  illustrate  the  application  of  a  new  definition  of  surface  density 
of  contamination,  the  results  also  illustrate  the  magnitude  of  the 
variation  in  exposure -dose  rate  to  be  expected  from  fractionation,  and 
the  sensitivity  of  the  results  to  several  important  variables.  During 
and  subsequent  to  the  preparation  of  thiB  report,  more  up-to-date  and. 
realistic  input  data  was  and  is  becoming  available.  This  data  includes 
total  chain  and  independent  fission  yields  for  various  fission  proces¬ 
ses,  improved  dose-rate  conversion  factors  from  new  decay-scheme  infor¬ 
mation,  and  fractionation-correlation  parameters  from  silicate -surface 
bursts.  More  detailed  calculations  will  be  mode  when  the  acquisition 
of  this  input  and  a  computer  program  for  its  application  have  been 
completed. 

For  the  present,  the  predicted  dose -rate  from  local  fallout  for  a 
given  degree  of  fractionation  does  not  appear  sensitive  to  the  inde¬ 
pendent  fission-yield  data  used  as  input.  The  sensitivity  to  total 
yield  is  significant  only  at  extreme  degrees  of  fractionation.  The 
predicted  effect  of  fractionation  varies  considerably  with  the  time  at 
which  exposure-dose  rate  is  estimated. 

Many  other  fractionation  effects  remain  to  be  dealt  with.  These 
may  be  of  much  greater  magnitude  and  show  much  greater  sensitivity  to 
independent  fission-yield  data  and  solidification  or  condensation  time 
than  does  the  effect  on  exposure-dose  rate  from  local  fallout  In  land- 
surface  burst  debris.  These  include:  the  partition  of  radionuclides 
between  local,  intermediate  and  world -wide  fallout  (as  treated  by  a 
theoretical  rather  than  a  semi -empirical  method);  the  effect  on  transi¬ 
ent  dose  and  dose-rote;  the  exposure  dose-rate  from  debris  which  is 
enriched  in  volatilely -behaving  mass  chains;  fractionation  in  venting 
underground  and  underwater  bursts,  on  land  surfaces  or  particle  sur¬ 
faces  that  show  different  correlation  properties  than  those  used  here, 
or  in  tower  and  low  air  bursts. 


REFERENCES 


1.  E.  C.  Freiling.  Fractionation  I.  High  Yield  Surface  Burst  Corre¬ 
lations.  U.  S.  Naval  Radiological  Defense  laboratory  Technical 
Report.  USNRDL-TR-385,  29  October  1959*  Also,  Science  133:1991 
(1961). 

2.  E.  C.  Freiling,  S.  C.  Rainey.  Fractionation  II.  On  Defining  the 
Surface  Density  of  Contamination.  U.  S.  Naval  Radiological  Defense 
laboratory  Technical  Report,  USNRDL-TR-631,  13  March  1963. 

3.  E.  C.  Freiling.  Fractionation  III.  Estimation  of  Degree  of  Frac¬ 
tionation  and  Radionuclide  Partition  for  Nuclear  Debris.  U.  S. 
Naval  Radiological  Defense  laboratory  Technical  Report,  USNRDL- 
TH-680,  12  September  1963.  An  abridged  version  appears  in  Science 
122:1058  (1963). 

4.  C.  F.  Miller.  A  Theory  of  Formation  of  Fallout  From  land  Surface 
Nuclear  Detonations  and  Decay  of  the  Fission  Products.  U.  S.  Naval 
Radiological  Defense  Laboratory  Technical  Report,  USNRDL-TR-425, 

27  May  i960. 

5.  E.  C.  Freiling.  Fractionation  in  Surface  Bursts.  Paper  in:  Radio¬ 
active  Fallout  From  Nuclear  Weapons  Tests.  Atomic  Energy  Ccmmis- 
aion,  TID-7632,  February  1962. 

6.  C.  F.  Miller.  Fallout  and  Radiological  Countermeasures.  Stanford 
Research  Institute,  January  1963. 

7.  C.  F.  Miller,  P.  Loeb.  Ionization  Rate  and  Photon  Pulse  Decay  of 
Fission  Products  From  the  Slow-Neutron  Fission  of  U235.  U.  S. 

Naval  Radiological  Defense  laboratory  Technical  Report,  USNRDL-TR- 
247,  4  August  1958. 

8.  R.  C.  Bolles,  N.  E.  Ballou.  Calculated  Activities  and  Abundances 
of  U235  Fission  Products.  U.  S.  Naval  Radiological  Defense  labora¬ 
tory  Research  and  Developnent  Report,  USNRDL-456,  30  August  1956. 


23 


9.  R.  D.  Present.  Phys .  Rev.  J2:7  (1947). 

10.  L.  E.  Olendenin,  et  al.  National  Nuclear  Energy  Series,  Div.  4,  9. 
McGraw  Hill,  1951*  Also,  L.  E.  Olendenin,  Massachusetts  Institute 
of  Technology,  Tech.  Rept.  No.  35  (1949). 

11.  K.  Way,  E.  P.  Wigner.  Piffil-  S£V.  ;£:13l8  (1948). 

12.  The  Effects  of  Nuclear  Weapons.  U.  S.  Atomic  Energy  Commission, 
revised  edition,  April  1962. 


\ 


APPENDIX  A 


COMPARISON  OF  OUR  NOTATION  WIB1  THAT  OF  OTHER  WORKERS 


The  following  comparison  of  notation  will  help  the  reader  compare 
work  published  by  different  investigators  on  the  subject  of  radionuclide 
fractionation  in  nuclear  debris. 

The  symbols  used  to  describe  fractionation  in  this  series  of  reports 
are  as  follows : 

F  Total  number  of  fissions  occurring  in  a  nuclear  event 

A^  Total  number  of  atoms  of  chain  i,  either  fission  product 

or  induced  activity,  occurring  in  the  event 

Yp  Total  yield  of  chain  i  in  the  event.  =  A^/F 

Number  of  atoms  of  chain  i  in  a  given  sample  as  determined 
by  radiochemical  analysis  and  corrected  to  the  time  the 
event  occurred.  The  same  symbol  has  also  been  used  for  a 
correlation  parameter  (see  b^  below). 

f^  The  number  of  fissions  in  the  sample  according  to  the 

analysis  for  chain  i,  hence,  the  mass-i  chain  equivalent 
fissions.  =  a^/Y^  =  F/A^ 

r,  ,  The  fractionation  coefficient  for  the  i  and  J  chains. 

‘->3  r  =  f  /f 

ri,j  V  J 

rog  ne  The  fractionation  ratio.  A  measure  of  the  departure  of  a 
given  samples  *  radiochemical  composition  from  representa- 
tivity. 

$  The  degree  of  fractionation,  <t>  =  log^rg^  ^ 


2b 


i 


b  The  slope  of  the  regression  curve  correlating  log  r, 

1  with  -0.  1  -  b^  is  the  slope  for  correlating  log  rp'gj 

with  *.  log  =  a1  +  (l  -  b±)<i > 

F  (t)  The  fraction  of  atoms  in  a  given  mass  chain  which  are  in  a 

R  refractory  (condensible)  form  at  the  time  t  under  the  fire¬ 

ball  or  cloud  conditions  prevailing  at  that  time. 

Edvarson,  Low  and  Sisefsky  (Al)  use  f^ys  which  they  call  a  "frac¬ 
tionation  factor,"  to  indicate  the  fractionation  of  nuclide  A  from  the 
mass-95  chain.  This  would  be  equivalent  to  our  rA,95>  if  we  1166,1  A, 
rather  than  i,  to  indicate  the  mass  chain. 

Mamuro,  et  al.,  (A2)  also  call  this  same  quantity  a  fractionation 
factor,  but  give  it  the  symbol  f. 

Miller's  notation  (A3)  is  more  complex.  He  uses  a  variety  of 
symbole,  which  he  calls  "fractionation  numbers,"  in  apparently  the 
following  way.  On  page  51  of  his  report,  he  uses  r0(A)  in  a  manner 
equivalent  to  r/v  q<5,  although  he  defines  it  on  page  xxvlii  more  like 
our  Fr.  The  quantity  RAy  on  page  62  of  his  report,  when  referred  to 
thermal  neutron  floulune  from  l£35  as  a  reference  event,  is  the  fami¬ 
liar  r99(a)  value,  frequently  used  by  radiochemists  and  diagnosticians. 
When  referred  to  a  given  nuclear  burst  as  the  reference  event,  however, 
it  becomes  r«  yy  In  our  notution  or  fA_ng  in  Edvarson' s.  On  pages  107 
and  IOQ,  rQ  (A, t)  la  our  FR(t),  while  r99(A)  is  again  rA;yy. 
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